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ABSTRACT ^ 

The equations defining errors generated during the spin-up of sym- 
metric rigid bodies and subsequent thrusting of the spinning body are 
derived from Euler’s dynamical equations and transformed into inertial 
coordinates. Examples of both spin-up and thrusting errors are given 
for the Ranger lunar landing capsule and an Apollo-shaped planetary 
entry capsule. 



The stability of nonrigid spinning symmetrical bodies is also dis- 
cussed. Equations describing the precessional motion are given, and a 
simple mathematical model is analyzed. Results indicate that to mini- 
mize errors in the spin-up and thrusting phase, a high spin-rate is de- 
sired, whereas to minimize errors during coast, a low spin-rate is 
needed. 




m 


I. INTRODUCTION 


More than a dozen space projects have successfully 
utilized spin stability as a mode of attitude control; many 
other space applications have been viewed analytically. 
This Report deals with three phases of spin stability — 
the spin-up process itself, thrusting the rigid body after 
spin-up, and the stability of not-so-rigid bodies during 
coast. 

Two projects presently under the auspices of the Jet 
Propulsion Laboratory (JPL) have investigated this area. 
They are the Ranger Blocks II and V and the Mariner 
Mars ’66. The Ranger rough-landing capsule was to have 
been separated from the parent spacecraft by spinning- 
up with canted nozzles (to a nominal 350 rpm in 1 sec), 
so that a positive separation velocity was attained during 
spin-up. Shortly after separation, a solid propellant 
rocket motor (which constituted roughly % of the total 


mass) was ignited, removing approximately 8800 ft/sec 
from the incoming velocity vector. It was necessary to 
determine the velocity vector at motor burnout. 

The possibility of a planetary entry capsule to be 
separated from the fly-by spacecraft was considered for 
the Mariner Mars ’66 mission, and is definitely being 
considered for future planetary missions. The physical 
separation of the capsule from the spacecraft, and the 
accuracy of the required velocity vector were studied. 
One method involved releasing the capsule from the bus 
with springs, spinning-up immediately, coasting for a 
while (to ensure minimal impingement of the rocket 
exhaust on the spacecraft), then igniting a rocket motor 
to give the required separation velocity. For most veloci- 
ties considered, a constant capsule mass (during burning) 
could be assumed. The question arose as to how accurate 


1 
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th* velocity vector was when added to the capsule, 
considering the spin-up as well as the thrusting errors. 
Another question which arose was with regard to the 
stability of the capsule after the entire separation ma- 
neuver (if it were still spinning at its design rate). If 
stability was required (say, for telecommunications), 
then a recurrence of what happened on Explorer 1 1 could 
not be tolerated. 

It is believed that the tools necessary for the solution 
of the three aforementioned problems, and sufficient 

Shortly after launch, Explorer L a satellite shaped like a cigar and 
spun about its axis of symmetry, began to increase its wobble angle 
until finally all spin was about the transverse axis. 


examples, are given so that the conclusion and equations 
may be generally applied. 

In the spin-up and thrusting Sections (III and IV), the 
applied torques are assumed body-fixed and constant. 
Neither gravity nor solar pressure affects the errors since 
they are inertial forces. Aerodynamic forces are not con- 
sidered. 

Since the areas discussed in Sections III, IV and V 
were investigated separately, each Section is self con- 
tained, and depends at most on the general discussion in 
Section II. 


II. RIGID BODY EQUATIONS OF MOTION 


A. Euler’s Dynamical Equations 

Consider the rotational analog of Newton’s second law 
of motion in inertial coordinates: 

t =s"> = l< nfl > “> 

where L is the external torque, / the angular momentum, 
II the inertia tensor 2 , and Si. the angular velocity. Re- 
calling that the time rate of change of a vector (A) in a 
moving reference frame of angular velocity (Si) with 
respect to some inertial frame is [^]i n «rti«i — [4 ] Mr 
+ [Si X A], Eq. (1) may be rewritten 

1 =[>*>],. = [> D, k +flx(nn) 

(2) 


2 For the results discussed in this Report, H is diagonal, i.e,, the 
cross products of inertia are zero. When a body of symmetry is 
designed to be spun, dynamic balancing can move the principal 
axes of rotation coincident with the geometrical axes to a degree 
beyond that of interest to this study. However, for the case where 
the products of inertia cannot be assumed zero, and for the effects 

of products of inertia on spin-up, an excellent treatment is given 

in Ref. 1. 


The SI in the above equations is that of the moving 
frame with respect to th° inertial frame. The torque (L) 
has been constrained to he body-fixed. 

Equation (2), when expanded, yields 3 ; 

I* fi) Wy <*)» — Lx (da) 

ly toy — (L — l x ) U>£ 0 >g = Ly (3b) 

/v to x (Zj ly) to x toy ” Lg '( 3 CJ 

Equations (3) are Euler's dynamical equations, the basic 

tool for the analysis carried out in this Report. 

If the z-axis is taken to be the axis of symmetry and 
of spin 4 (whereupon l x ~ l y ), then a convenient method 
of presenting the cross-angular velocities is obtained by 


3 The rate of change of inertia component in -^L (II £2),i.e., 

etc., must be assumed zero, since L, etc., must necessarily be 
caused by a mass expenditure which is no longer part of the rigid 
body. If mt, and L, etc., is not the result of a dm/dt , then the 
assumption of a rigid body is no longer valid. 

4 Spin is defh ed as the rotation about the body roll axis (<*>,) only, 
and the cross plane rotation («, and w*) is referred to as wobble. 
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letting o) = w, + i o> yj thereby mapping both components 
onto the complex plane. Noting that w — w* + i w y , Eqs. 
(3a and 3b) combine to give 

o) ~~ i A o) w* == N (4) 

where A - UL - 1, (- 1 < A < 1), and N = h±lhi 

*X 

(cross-angular acceleration). . ui^e l £ = I y , Eq. (3c) re- 
duces to 

i) z = (angular acceleration of spin) (5) 

It will be shown later that N~ must be a constant for any 
reasonable closed form of Eq. (4) and subsequent equa- 
tions. 

If all torques are zero or the body, then <d z — con- 
stant = So (initial spin rate), whereupon Eq. (4) is 
immediately integrated to 



LINE OF 
NODES 


w — (o 0 C 


( 6 ) 


Fig. 1 . Euler angle transformation 


where a> 0 = w 0 x + i o> ( , u . The components of Eq. (6) are 

Wj — w woj COS A So t <i> 0 {/ sin A Sq t 

0)y 1 3 W Wo y COS A Sq t ”f WQjp Sin A $0 ^ 

When time is eliminated as a parameter, it is seen that 
the magnitude of the cross-angular velocity i ' constant, 
w = (%ir w ow) ,/<J > an< ^ routing about the angular mo- 
mentum sector / at a frequency A s 0 . 


Note that the magnitude of the spin rate, w*, is de- 
pendent on s 0 and N z alone, hence would remain con- 
stant regardless of f he value of N = N x + i N v . This is 
true only because the body is symmetrical. 


shown in order on Fig. 1, and a^e discussed more fully 
in Ref. 2. They are 




cos <t> sin ij) 0" 

— sin <f> cos <p 0 ; 

0 0 1 




"1 0 0 " 

0 cos 0 sin 0 ; 
_0 —sin 0 cos 


w 


cos \j/ sin 0“ 

— sin ^ co s^ 0 

0 0 1 _ 


( 7 ) 


B. Inertial Transformation 

The transformation from body-fixed coordinates to 
inertial coordinates is made with conventional Euler 
angles.’ 1 The X 0 — To Z# frame is inertial, the x - y — z 
frame being body-fixed. The three angular rotations are 


As seen in Fig. 1, the transformation from Euler rates 
to body-fixed rates is 


w* 


"sin 0 sin ^ cos \f> 0“ 



Wy 

= 

sin 0 cost)/ —sir up 0 


0 



cos 0 0 1_ 


J. 


’Anothei Eulerian angular transformation for analog computer sim- 
ulation is given in Appendix A. 


The full angular transformation [\p] • [o] • [<*>] is discussed 
in Section IV-E. 


3 
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If the cross-angular \ elocities (w r -r i a v ) are mapped 
onto the complex plane from Eq. (8), 

iu = (0 -f i $ sin 0) e ,e (8a; 

which becomes (when 4> is eliminated): 

w = [0 -4- i (w z — <f) tan 0] e ^ 

If a small angle approximation for 0 is made, which is 
valid in light of the small errors to be encountered, 
then (<? «=; tan 0) 

G> = £0 *r i (w r — \j/) 0j e~ 

Defining a complex angle of attack by 
a = Be ^ 

then » and a are related by 

a + i o» r a = w (9) 

which can be verified by substitution. This complex angle 
of attack a — 0 (cos * — i sin \5r) is seen from Fig. 1 to be 
the angle of attack (0) as seen by an observer sitting on 
the X 0 -axis for a and the — Y 0 axis for 3m a. How- 


ever, a is only a mathematical intermediary to obtain tTe 
inertial angle of attack in complex form. This is given by 

a, = a (l r ? 

This cr.n be shown as follows: From Eq. (8), <u c = 0 cos 6 
+ ^ and cos 0 ^ 1, 

i J u>; dt = i cos 0 + ^) = i f (cos 0 

= i (<f> 4- i/r) 

whence Eq. (TO) becomes 

a, = (XI) 

Equation (31) then is the full inertial transformation, for 
expanding fl, = 0 (cos ^ + i sin 4>), which is the angle 
of attack 0 as seen by an observer sitting on the X 0 -axis 
when <f> = 0, and the Y 0 -axis when <f> — r/2. Thus, SRf a t 
will give the X 0 -eomponent of the time-varying angle of 
attack, and 3m dt, the Y 0 -component. 

In a similar manner, it can be shown that the complex 
angular rate as measured in the inertial X 0 — Y 0 direc- 
tion is 


III. SPINNING-UP OF SYMMETRIC RIGID BODIES 


When a symmetric rigid body is free in space with an 
initial rotation (0 O ) about its pitch-yaw axis, and is then 
torqued about its roll axis to a spinning condition, three 
errors result: 

1. The angular displacement at the moment of spin-up 
due to the initial conditions 0 O and 0 O . That is, with 
respect to the inertial reference, there is an angular 
displacement error (0 O )> a rate error (0 O ), and a time 

from release to spin-up (8f) which causes an angle 

of 0 O + 0 o 8t between the body-fixed z-axis and 

the true inertial axes at the time between release 
and spin-up of the body. Since 0 O and 8f are 

parameters which are totally independent of this 

analysis, their effect is omitted here, but must be 


determined and vectoriallv summed with the errors 
listed below for a full separation error analysis. 1 * * * * 6 
The new inertial axes are then defined as the body- 
fixed axes at t = 0, thus omitting this angular error 
from further consideration. 

2. When the body is spun-up about the new set of 
inertial axes, it is found that the angular momen- 
tum vector, which defines the direction of spin, is 
displaced at an angle a,. This rotation of J is the 
result of both the initial tumble rate (0 O ) and the 
errors associated with the spin-up itself, which arise 
out of the body trying to spin-up about axes other 


6 See section III-D. 
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than principal, thus generating cross products of 
inertia. 

3. The body-fixed z-(spin) axis rotates about the 
angular momentum vector at an angle p. This wob- 
ble is caused also by both initial and spin-up errors. 
In the analysis, p is derived first since the develop- 
ment of 0 / depends on p . 

A . Spin Jet Misalignment 

The errors in the spin-up are caused by an initial 
condition and an error associated with the alignment of 
the body-fixed torque vector. This error in the torque 
vector is made up of the five components shown in Fig. 
2. If there are n spin jets, the nominal thrust of each 
being F, and if the jets are evenly spaced about a circle 
of radius R, the plane of which is i from the c.m. of the 
body, thr torque in the z-direetion is (using small angle 
approximations) 

L z - nFfl (12) 

The possible effects of AR on L z are neglected. By in- 
spection of Fig. 2, the RSS cross-plane (i.e., x-y) torques 
are 7 

L = jnP [(my- + (My + (,l) s ] + n(MF)=}* 

(13) 

If AF = kF, which is the case when spin jets are mani- 
folded L = yfn F (a -4- ib), the complex number being 
for the direction of the net cross-plane torque which is 
not known from Eq. (13). However, when the cross- 
plan e torque is compared to the primary torque, it is seen 
that 

L — L x + i L v = y L z (14) 

7 Assuming the misalignments are randomly distributed and no 
biases exist, especially in alignment tooling. If biases do exist, then 
Eq. (13) is conservative. 



Fig. 2. Spin jet misalignments: a. top view; b. side view 


where 

a *f i b 
7 ~ V'nfi 

Equation (14) is the form which wall be used for the 
cross-plane torque, y is defined as a real number without 
loss in generality. Furthermore, the body-fixed x-t/ axes 
and the inertial X 0 -Y.i axes are defined to make the ini- 
tial rate 6 it and thf cross-plane torque L act in the same 
direction, thus pr renting a worst case. 

B. Primary Torque Considerations 

The cross-plane angular acceleration is 

N = T = = T yN -~ = < 1 + < 15 > 

ij- It 

Eliminating from Eqs. (4) and (5), 

= y(l + X)N Z (16) 

It is immediately seen that the form of N z governs the 
solubility of Eq. (16) 8 . In general, the / N z dt utroduces 
a constant of integration, the initial spin rate m oz . If par- 
tial restraint of the body during spin-up is made, <u 0 z > 0, 
and Eq. (16) has no analytic solution. Therefore, w O r = 0 
from here on. 

Even after making = 0, the form of N~ is quite 
restricted for closed solutions. If a cold gas spin-up sys- 
tem is considered, then j N z dt — c( 1 — e~ H ), since an 
exponential decay of thrust (and an exponential increase 
in spin rate) is the best approximation, Av can be verified 
by substitution, the constant value c acts in the same way 
that wo* does in the solubility of Eq. (16), preventing a 
closed form solution. Polynomial approximations to the 
exponential either introduce constant terms like the one 
above, or consist of functions which have singularities 
at t = 0. However, since the exponential thrust spin-up 
has a maximum torque early in the spin-up (decaying to 
zero) the errors associated with this form are less than 
those of constant torque, and the method presented here 
is conservative. Valid approximations to the exponential 
with the constant thrust method can be made. Hence, 
it may be concluded through some devious logic that 
N z must be constant for a closed form solution of Eq, 
(16) and subsequent differential equations, and 
ui S = J* N z dt — N z t . 

8 If a digital solution of Eq. ( 16) and subsequent differential equa- 
tions is used, any representation of 3V S is allowable. The statements 
made here pertain to more ‘‘closed” solutions, in terms of known 
functions, which also require digital solutions, but can be easily 
generalized. 
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C. Wobble Error During Spin-Up 

The angle at which the body-fixed roll axis wobbles 
about the angular momentum vector can now be deter- 
mined. When N z is constant and w 0 * is zero, then Eq. (16) 
is integrated to 


(spin), 11 and J. When t = t x and all external torques have 
ceased, J is in valiant in space while the spin-axis * and the 
angular velocity vector 11 coplanar with J and % rotate 
about J at an angular rate XN.fJ. The angle fi, which is 
the wobble angle between * and J is seen to be 


u.v, 

- y N z (1 + a) e~ 


i: 


t* j , dd 
e 2 dr -f cu 0 e 2 


(17) 


where cd 0 = <u 0 * 4 i <a ov , and 0 < t < t l (torque time) is 
the running variable. The integral in Eq. (17) is a Fresnel 
integral, the results of which have been tabulated after 
the proper change of variable (Refs. 3 through 7). 


tan | p j 


VTT+J ; _ /, |<» 

J\ I' So 


or, letting tan j p ! « j f) j. 


1*1- 


(1 + X) So 


(18) 


The angular velocity vector, 11, may then be written 
in body coordinates 


11 — w* £ + ci>y j -f , 


where w- = 3l< <*>, w y = 3 m « from Eq. (17), and w* = N z t. 
When* = t x (the upper limit on torque time), o> will assume 
the value attached to Eq. (17) when fj is substituted for t , 
and w* ~ N s t x = s 0 (the final spin rate), a design quantity 
which is sure to be known. The angular momentum vector 
is then J = 1 1 1, where I is the column matrix F = ( I x I y 1 2 ). 
Figure 3 shows the relationship among the body z-axis 


In complex form, Eq. (18) is 


(1 4* A) So * 

which is seen to be upon expansion using Eq. (17) (£ = * 3 ). 


p = ^ e™* 
S 0 


'f 


-i A.V 


-T» 7 , 

t? 2 ttr + 


Letting 


and 


(1+A)$ 0 

N; iHl ,;/■'* zdd. T . . 

= — « 2 I e 2 dr 


e 2 1 


K, = — - e * 

So 


(19) 


z 



Fig. 3. Wobble angle, /? 


1 

= (1+A)s 0 e 2 (20) 

then 

/3 = y Kj + <i> 0 K 2 (21) 

where K t and K 2 are functions of the shape of the rigid 
body (a) and of specific spin-up values (N t and f,), but are 
independent of the error producing parameters y and u> 0 . 

Looking at Eq. (18) again, 



y N, f** 

- — / e 2 
«o Jo 


T ’dr + 


top 

(1 4* A) $q 



e 2 


dr + 


to o 

(l + A) So 


( 22 ) 


Equation (22) is the form of the wobble angle which 
shall be used. The ] /J | given in Eq. (22) is the initial 
condition for the thrusting phase (to be discussed in Sec- 
tion IV). 
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D. Inertial Displacement of J During Spin-Up 

The angle through which J moves with respect to the 
inertial X 0 -Yo-Z<> axes is now determined. Recall that Z 0 
and z are coincident at t — 0. When torquing begins, 
J builds up from 0, wobbles about Z„ then when t = t u 
J comes to rest at angle a, from Z 0 . Consequently, when 
determining a t , the inertial transformation defined in 
Section II is needed. Using the value of u found in Eq. (17) 
in the inertial Eq. (9), 

. «a - v, „ ft -iky t , 

a 4- i N; t a — y N x (1 +A) e 2 I e - T dr 

the solution of which is 

a e 2 * = a 0 + K 3 y 4- K 4 u> 0 (23) 

where 


after spin-up (e.g., thrusting), then a t and p must be 
combined in either of two ways: 

1. If there is a coast period between the end of spin-up 
and start of maneuver, so that the position of * with 
respect to / may be considered random (i.e., % has 
precessed many times about /), then a t and fi are 
independent and treated as separate errors; i.e., a f is 
the error incurred during spin-up and p is an initial 
condition imposed upon the subsequent maneuver. 

!«r=M*+ !/?!* ( 27 ) 

2. If the maneuver takes place immediately after spin- 
up, such that the position of * will be known with 
respect to J, then a t and p must be vectorially added. 
In this case, the spin-up parameters must be well 
defined, and experimentation must agree with theory 
so that the directions of a t and P in fact correspond 
to the analysis. 


K 3 - 


/ t i <1+A)JT, ( fT -l‘A.V y ) 

e~T~ T2 j j e ~T~ c dddr 


and 


K 4 = 


ft iS- (i+A) , 

= i e 2 


(24) 

(25) 


The integration in K 4 is seen to be of the same form as 
the integral in (Eq. 19), i.e., a Fresnel integral. 

Then, from Eq, (10) 

ij* W 3 dt ifs r tdt 

at = ae = ae = ae 


and Eq. (23) becomes 


fl/ = Oo + K 3 y + K| <i) 0 (26) 

Since at t = 0, z and Z 0 are coincident, Oo = 0. However, 
for an overall analysis, it must be remembered that 
Oo = ©o + Oo 8f (discussed in the Introduction). 


E» Combination of Body-Fixed and Inertial Errors 

If the spin-up error is to be considered by itself, i.e., not 
used as an initial condition for some maneuver after spin- 
up, then the error in the angular momentum vector a 7 is 
considered the prime error (p being a wobble super- 
imposed on a,). However, if a maneuver is to be performed 


It will be shown later that the wobble rate <£^(1 + x)s 0 , 
so that as A — > — 1 (a long thin rod spun about its axis 
of symmetry), <f> can become small, saying that % rotates 
very slowly about /, hence the direction might well be 
predicted. For the vectorial addition, refer to Fig. 4. 


z 0 



Fig. 4. Combination of a 7 and p 
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Since a t = 0 (cos <f> 4 - i sin <J>) when <f> = 0 (standing on the 
X 0 axis), an observer sees iR* «/ — 0, consequently it is just 
that component which is then projected onto the re- 
direction (in this case — 0, since 0 is defined counterclock- 
wise). The same is true for J) m a,, but no sign change is 
required. p is defined in the normal way. Hence, in 
Fig. 4, a unit sphere is superimposed on the vector dia- 
gram, and the angular components are shown projected 
on the sphere. It can be seen that 

€ = P fRr Ck i 


which can be verified by substitution. The final form of 
the K’s is reached by getting rid of the clumsy design 
parameter N : . Note that, for constant N~, 

N 2 = sjt x 

Then, by substitution 

and the four K s assume their final form 


Sie e — 3lep + 3mdl 


whence 


€ = p id/ 


(28) 


and j € j is the required error. 


A discussion of the errors defined above follows, with 
examples given for e, p and a t . 


F. Discussion of Integrals 

In order to solve numerically Eqs. (17) and (26), the 
quantities K u K 2 , K 3 and K 4 must be evaluated. In all 
cases the variable is t where 0 < t < t x . As mentioned 
previously, K x and K 4 are Fresnel integrals, and the double 
integral in K 3 is the integral of a Fresnel integral. To get 
these integrals in familiar form, a change of variable is 
required. The usual argument in the exponent of the 
Fresnel integral is brtr/2. Hence, if this substitution is 
made in K 3 and K,, they become 



r*r XT' 

r \ * 

I Z — fl / . . 

= M • \ e ~' U 


K,(t) = 


(1 + A) s 0 


iX Jo ft 

e 2 


(29) 

(30) 


, WS- 

K ,(*0 = T I e ^“ 2 / e dxdu 

J 0 ,/o 


«.<*•) = 


V #o ( 1 +A) 

— ~ 

I / — « 2 /h, 

\ (1 + A)s„ jf S 2 


(32) 


When the shape of the body is stipulated, and the spin-up 
parameter s 0 is given, then the values of K = f (t x ) may 
be determined, and the errors may be calculated with 
assumed values of y and <i> 0 . As an illustration, Figs. 5 
A through 14 show the K’s as functions of t x for a few com- 
binations of A and S 0 . 




The Fresnel integrals of K, and K t are of the following 
form (see, for example. Ref. 4): 

C (o) = J cos — u a du 
S (o) = y sin y u : du 
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where 


E 



e 2 


du = C (k) 4 


iS (u) 


may be formed, and the conjugate of E is 


du = C (c) - iS (v) 

It is seen immediately that C(0) = S(0! = 0; however, 
C(v) and S(c) must be expanded asymptotically to de- 
termine 



lim C ( o ) = lim S (v) — -i- 

t +'jc V — » %> “ 


Reference 3 discusses briefly the nature and expansions 
of the Fresnel form in K 3 . However, no tables of general 
parameters are known to exist. When the above proper- 
ties of the Fresnel integral are applied to the K’s, it may 
be deduced that 


K.WKO.95^-^ 

(33) 


(34) 

K 3 (ti) | ^ k*iT 

(35) 

| K ,(«,)|< 0.95^ (1+xs 

(36) 


where the constants are the maximum absolute values 
of the oscillating terms. 


G. General Observations 

The K’s are bounded from Eqs. (33-36); therefore, the 
following conclusions (some obvious) may be drawn 
about the errors. 

1. As t x increases, keeping s () constant, 


♦Where k is some finite number. This can be seen from Eq, (31), 
where the inner integral is a Fresnel integral, which is always 
bounded. This Fresnel integral is multiplied by a bounded sinusoid 
c iiru */2 anc ] integrated, the result of which must be bounded. 


t j K x (t x ) j decreases, implying that the error con- 
tribution of the effective spin jet misalignment (y) 
to the wobble angle (ft) deer ases as spin time t x 
increases. 

ii I K 4 (f t ) i increases without bound, stating that 
the initial tip-off rate, w„, has had an increased time 
to act through, thus increasing a h 

2. As the value of s 0 is increased, keeping t x constant, 
K u and K 4 decrease (K s not being affected) This 
indicates that s 0 should be as large as possible, 
which could be expected from the gyrodynamic 
effect of increased spin rate. 

3. The quantity ] K 3 (t x ) j is independen* of t x and s 0 
in its limiting condition (Eq. 35). K 3 contributes to 
the a l error as amplified by y. A physical interpre- 
tation of this is not immediately apparent. 

4. All four K’s are functions of the shape of the rigid 
body (A). 

5. As t x -* 0, the four K’s reduce to (using L’ Hos- 
pital’s rule on K x ): 


K x ( 0 ) = 1 


K 2 (0) = 


1 

(1 + A) s 0 


K 3 (0) = 0 


K, (0) = 0 


which when substituted into the equations for a h ft 
and then e yield (ai = 0), 


€ 


0 = Y + 


(1 “ A) S 0 


Jz + 

h 


(37) 


which is immediately seen to be the rotation of J 
caused by the instantaneous application of the spin- 
jet misalignment, plus the initial tip-off component 

L/L. 


Equation (37) can be used as a quick method of approxi- 
mating total error. However, it should be remembered 
that this idealized case (i.e., t x — 0) is not conservative at 
all, and usually represents the lowest value attained. 
However, if t x is small, and an appropriate fudge factor is 
used, Eq. (37) is an * order of magnitude” approximation. 


9 



JPL TECHNICAL REPORT NO. 32-644. 


H. Limiting Cases 

1. A = 0 (sphere). The values for the K s are (using 
L'Hospitals rule on K x and K) 


Kj = 1 

“ 1/ #0 



and the errors are 


O ~ -L W ° 

0 ~ 7 + 

So 



€ = ft — 10/ = €* + t€ y 



Notice that the \/^ in the numerator of the second 
term of a 7 makes e divergent when <o 0 0. As t, 0, 

€ = /?, the value of Eq. (37). 


2. A = — 1 (thin rod). As A — > — 1, K 2 increases without 
limit (fC 4 is finite because (1 + A) appears in the upper 
limit of the integral), the other Ks being bound. 
This means that spinning-up a long thin rod about 
its symmetry axis (l z = 0) when an initial tip-off 
rate is present, is highly unstable. This is intuitive 
in that if l z — 0, there is no spin momentum to 
counteract the tumble of the rod. 

3. A = + 1 (flat disc). As A + 1 (its maximum value), 
the coefficients of the Ks are minimized, indicating 
that for a given spin-up time and rate, an infinitely 
flat disc is optimum. This is obvious from gyro- 
dynamic considerations. 

I. Initial Conditions 

t The relationship between the initial tip-off rate (w 0 and 
d 0 ) is found from Euler’s rate transformations (Eq. 8) 

• * 
u> # sin 0 sin i p + 0 cos 


» • 

o:fj — <f> sin 6 cos \f/ *- $ sin \p 
• « 

<D Z ~ <f> COS 0 + if/ 

Referring to Fig. 1, it may be stipulated that the initial 
conditions at t = 0 are 


= 0 

(1 

o 

<9 = 0 

i: 

\p *= 0 

O 

11 


which result in 

toor — do J 0 j — 0 

where 

<i->UX ”f* 1 O) Qy d 0 . 

Thus, it is being stated, without loss of generality, that 
whatever direction the tip-off occurs in, defines the X o -Y 0 
directions. 

J, Numerical Evaluation of Errors 

As mentioned in the Introduction, the two cases of 
interest in this Report are the Apollo-shaped capsule and 
the Ranger landing capsule. Since the K’s need only the 
body shape A ana final spin rate s 0 to be evaluated, it 
was decided to keep the error sources d 0 and y strictly 
parametrical. Spin-up time (£ x ) is taken with wide enough 
bounds to encompass most situations. The values used in 
determining the Ks are the following: 



APOLLO 

RANGER 

A 

1/3 

-3/4 

Final spin rate, rad/sec ($,,) 

1, 2, 5, 10 

30 

Spin-up time, sec (t,) 

0-10 

0-3 


Plots of both the real and imaginary portions of Ki, K 2 , 
K and K 4 vs t x are given in Figs. 5-14. 

The numerical evaluation of K u K a and K 4 was, at first, 
very frustrating. Fresnel integrals are known for their 
difficulty in approximating over a large range of argu- 
ments. However, the excellent orthogonal polynomial 
(Chefcyshev) approximation given in Ref. 6 is satisfactory 
for K and K 4 , while the corresponding function-generating 
polynomials were of sufficient accuracy for K 3 . Asymp- 
totic expansions for Fresnel integrals and for integrals of 
Fresnel integrals are given in Ref. 4. 



Im (K) 




Fig. 8. Imaginary parts of K u K 2 , K and K 4 vs fi 
(s 0 = 2 rad/sec, A. = + Vi) 


For Ranger, the specification for y is 0.006 rad with a 
spin-up time of 1.0 sec. 0 If an initial tip-off rate of 
(o — 0 O = 0.005 rad/sec is assumed (realistic), then from 
Figs. 13 and 14, 


Fig. 6. Imaginary parts of K„ K :r K : , and K t vs fi 
(s 0 — 1 rad/soe, \ — + V») 


°Aeronutronic Publication U-200, Final Technical Report— Lunar 
Rough Landing Capsule Development Program, Newport Beach, 
Calif., pp. 3-28. 


It 










0 2.0 4.0 6 0 8.0 10,0 0 2 0 4 0 6 0 8 0 10.0 

t \ , sec /, , sec 


Fig. 10. Imaginary parts of K u K 2 , K 3 and K 4 vs fj Fig. 12. Imaginary ports of K u K : , K, and K, vs t, 

<s 0 = 5 rad/ sec, k = + V* ) (s„ = 10 rad/sec, A = + V*) 


(Re K, = - 0.208 3le K 2 = + 0.034 

(Re K, = - 0.238 (Re K 4 — + 0.263 

Am Kr—+ 0.094 An K 2 = + 0.129 

An K> - + 0.671 Am K 4 = + 0.438 


where 

(3 = — 0.0011 + 0.0012 i and a, = - 0.0001 + 0.0062 i 
and e = 0.005 + 0.001 i and | « | = 0.005 rad 
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'" MC (s„ = 30 rad/sec, A = - 3 A) 

Fig. 13. Real parts of K u K 2 , K , and K 4 vs t, 

(s„ = 30 rad/sec, A = — %) 


IV. THRUSTING OF SPINNING SYMMETRIC RIGID BODIES 


A great deal of work has gone into analyzing the effects 
of rocket thrust misalignment on accelerated spinning 
bodies. References 4 and 8 are the classic works, done 
during the 1940's, primarily on spin-stabilized rockets 
launched from launchers with and without fins. The 
basic theory used in this Section is a combination of these 
and other referenced works. No claim is made to origi- 
nality, except, perhaps Eq. 67. 

The one design parameter which is common between 
the last Section and this Section (except shape of the 
body, A) is the design spin rate s n , An output of Section 
III, the wobble angle /3 n > along with s {h is used as an 
initial condition here. The inertial axes X (l -Yo-Z rt are rede- 
fined for this Section. When t = 0, the Z 0 axis coincides 


with the angular momentum vector, J'°, and the Y 0 axis 
is along the direction of — J3„ (coinciding with the 0 direc- 
tion when <£ = 0; see Fig. 1). 

First to be determined in this Section is the angle of 
attack, a h of the thrusted rigid body 11 . This inertial angle 
of attack defines the position of the thrust vector as a 
function of time. The thrust vector is then integrated 


10 Thus, the error between Zn and J of Section III, a !y is subtracted 
out here. Of course, this error must be considered, and is dis- 
cussed in Section VI. 

n This ai differs from that in Section III in that (a), this a f is in a 
different time domain since the Xci,Yo,Zo axes are ledefined when 
£—0 and (b), the at here includes wobble angle. Hence, o/ () 
when f = 0. 
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into velocity, and the resultant error in both the magni- 
tude and direction is considered. There are certain sim- 
plifying assumptions made throughout this Section, such 
as constant mass, inertia, etc. A discussion of the effects 
of these assumptions is found in Appendix B. 


here, and p, a measure of the misalignment, is estimated 
by £ to the accuracy of small angle approximations. 
Appendix C contains data on the effective thrust misalign- 
ment from the 2nd and 3 rd stages of the Scout launch 
vehicle. This data was used in the numerical evaluation. 


A . Thrust Misalignment 

Thrust misalignment, which can be defined as the 
lever arm through which a given thrust acts, gives the 
body both a rotational and a translational acceleration. 
The translational motion will be neglected in this Report 
since in nearly all applications of the dynamics described 
herein, such motion is second order. Therefore, the need 
for considering two separate types of thrust misalignment 
is avoided (see Fig. 15). If F is the thrust vector, constant 
with respect to the frame shown, 8 the linear displacement 
of the point of thrust application (e.g., due to uneven 
rocket nozzle erosion), and a the angular misalignment 
(due to rocket nozzle misalignment), then the bver arm 
(p) through which F acts is seen to be (letting a « sin a), 
S + ia. If, however, an “effective thrust misalignment” 
angle £ is defined such that £ = S/i 4- a, then the lever 
arm is 

P = u (38) 

which requires only the definition of the effective thrust 
misalignment angle, £. This assumption is allowable, since 
in the true definition of thrust misalignment, part of the 
a error went into translational motion which is neglected 



It is also required that the thrust vector intersect the 
roll axis, i.e., no roll torque is present, if this is not the 
case, then nonconstant effects on spin rate during burning 
must be considered, which make the differential equa- 
tions nonlinear (see Appendix B). 

In determining the torque vector from Fig. 15, 

F = t ?fc 


whence torque 


L = rXF = lF£i 

(40) 

which is all in the x-direction. As stated before, all the 
quantities constituting torque (i, F, £) are constant with 
time. 

B. Body-Fixed Equations of Motion 


Substituting the above torque into Euler’s 
equations (3) gives 

dynamical 

55 

1! 

1! 

3 

3 

! 

•3 

(41) 

w z y~ — 0 

(42) 

Letting / oi z dt = s 0> the solution of Eq. (41) is 


/ W \ iXf , , iN 

<0 = ( OiQ ) £ X V + — 

\ X So / X 

(43) 


where O>0 — (Oox + iu)o v is discussed below. 


The wobble angle, /?, as a function of time, is given by 
Eq. (18) as 

«*> 

Resolving Eq. (43) into components, noting that Sm N = 0, 
g)jj “ iRf (i) u>o x cos \s„t ” (woy N/ A®o) sin Ks 0 t 

otp m w (1)01 Sm \$ot 4 “ (u>oi/ " N/ A$o) COS As ot Ai / As 0 

which reduce to (when time is eliminated) 


Fig. 15. Thrust misalignment angles 


( O * + (top — N/KSq ) 2 — < Og s + (wok — N/ASo ) 2 — R s (<i>o) 

(48) 



II 
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Fig. 16. Body-fixed rates of a symmetric body 
with torques 


which is a circle with its center at (0, N/\s 0 ) ? and a radius 
ofR(a) 0 ) = 4- (w 0y — N/\s q ) 2 }±. The complex radius 

vector, o)(t) rotates around the circle at an angular rate 
A So (see Fig. 16). Since the coordinates are body-fixed, 
the angular momentum vector (J) and the angular ve- 
locity vector (£2), which are coplanar with 15*, rotate 
about the axis of symmetry of the body. 


C. Initial Conditions 

Certain relationships among the Eulerian angles need 
to be derived before initial conditions can be stipulated. 
Just before ignition of the rocket motor, the spinning 
body has no torques acting upon it, and the form in which 
the initial conditions appear must be determined. It is 
given that when t = Q, # = # 0 , <o oz = s 0 > 0 = S 0 = ft 
(from Section III), and <j> = <j> 0 (precession rate). Needed is 
wo = wox 4- iw oy in ternr, of the given quantities. Without 
loss in generality, the other Eulerian terms 0, <j>, $ are 
zero initially. The proofs that 0 = 0 and = constant v/ill 
not be given here. The Eulerian angular transformations, 
Eqs. (8) are 

o>or ~ 0 

* * 

woy = " <f > o sin So ^ <j > o 0 O (49) 

• • * • 

So — Wo s ~ <£o COS 0 O 4- \j/ 0 *=* </>o 4- \f / o 

whereupon 

<i>0 *” 4“ t <0 oy 1 <f > 0 $0 

Since it is the second order differential equations 
(<o = d 3 $/dt 8 , i some angular displacement) that are of 
interest, only two initial conditions must be stipulated 
($o, G u ). Hence, some relationship must exist between 
and So and Co. This can be found by eliminating <o in Eqs. 
(4) and (8a) when t = 0 and N = 0 and initial conditions 
exist. 


to — i A <o<o® — N --- 0 (4) 

<o = (0 •, i t^ sin 6) e-W (8a) 

Taking the time derivative of Eq. (bu) and inserting the 
initial conditions, i.e., 6 = 0 o ,<j> = = fa, 8 = 0, £ = 0, 

and \fi = 0, 

$ sin 0 U (<j/ 0 + A. to 0 ®) = 0 (50) 

which has roots of <£ — 0, $ 0 — 0 and = -A to 0 ®. If 0 O = 0, 
there are no initial conditions (since there is no wobble). 
If 6 0 7^ 0, which necessarily implies that </> 0 0, then 

ij / o = — Ac.;,,.- for Eq. (50) to hold, and when this value for 
^o is substituted into Eq. (49), 

* _ ( I + A) S 0 
COS $0 

and Eq. ( 9) becomes 

Wo = <^0^0 = 4- A) So 00 

for small 6 0 . 

D. Complete Inertial Transformation 

Now the value for a ; (the inertial angle of attack) can 
be determined. Equation (9) stated that a 4- i s Q a = w, 
and when the value for (found in part B) is substituted 
and the differential equation is solved, the result is 


(51) 

(52) 


= “* + {.1XTT) + Ifjf+ij} { 1 " 

(53) 

and from Eq. (10), 

Ci — a = a e' iqt 


Letting <o 0 — t (1 + A) s 0 6 0 and a 0 = a, 0 ~ 0 o , Eq. (53) 
becomes 

- = «• + - *•} { l - 

(54) 


Without going into too much detail at this point of the 
development, it can be seen that a, is the sum of two 
sinusoids of differing amplitudes and frequencies. The 
frequency s 0 is the spin frequency, and (I + A) s 0 is seen 
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from Eq. (51) to be <£ 0 , or the initial precession rate if 
cos 0 P ^ 1. The amplitude of the first term can go to zero 
if N counteracts the wobble angle, 0... However, this only 
nulls the amplitude to the precession part, and does not 
cause 0 / to go to zero. If 0 O = 0, the amplitudes of the two 
sinusoidal terms differ by a factor of (A -f 1), which is 
bounded by 0 and 2. It will be seen later that if 0„ ^ 0, 
then it can become a predominate term. 

E. Inertia I Thrust Vector 

Since the inertial angle of attack («/) is now known 
with respect to body-fixed parameters, the motion of the 
thrust vector in inertial space can be determined. The 
tool needed is the full Eulerian angular transformation. 
Rewriting the rotation matrices (see Fig. 1), 


vector (F) is given in Eq. (39). If a small angle approxi- 
mation is made for 0, then the inertial thrust vector, F', 
resolved into its coordinates from Eq. (55) above, is 

F' = — F £ (cos <f> sin $ + cos ^ sin <£) + FO sin <j> 

F' — — F £ (sin </> sin $ ~ cos <f> cos $) — F 0 cos <f> 

F' = F£ 0cos # F cos 0 (57) 

The cos 0 is retained in F' for later consideration. 

If a complex inertial thrust vector is defined F' 

= F' + iF\ then bv substitution, it is found that 

F 'ry -fe*} (58) 

The exponent in the first term is determined from (cos #~1) 


M = 


M = 


M = 


[ 

[ 

[ 


cos <f> sin <f> <fj 

— sin <f> cos <f> 01; 

0 0 lj 

1 0 0 I 

0 cos# sin# I; 

0 —sin# cos#J 

cos ip sin ^ 0 

— sin ^ cos ^ 0 

0 0 1 . 


J i w : dt — i J (<j> cos # + ip) dt — i(<j> + iji) — s 0 t 

(59) 

and the second term is by definition (Eq. 11), a,. There- 
fore, 

F' n = » F {{ e'*"' -a,} (60) 

The above equation is the desired form of the complex 
inertial thrust vector. It is seen that if the two error 
sources { and # 0 are zero, = 0. 


Then the operation of [^] • [#] • [<£] = Ton any inertial 
vector will transform the same vector into body-fixed 
coordinates, or conversely 


= r* A 


body- fixed 


j>] • [#]■'* W~ X Shorty- fixed 

(55) 


where the matrix multiplication yields 


r -i = 


cos cos ^ — cos 0 sin <f> sin ^ 
sin ^ cos ^ + cos 0 sin ^ cos <f> 
sin 0 sin ^ 


— cos^sin^ — cos0cos^sin<£ 
—sin ^ sin <f> + cos 0 cos <j> cos ^ 
sin 0 cos ^ 


sin <f> sin 0 " 
- cos <f> sin 0 
cos 0 


(56) 


Since r~ l is orthogonal, i.e M j | = 1, the inverse of the 
matrix is equal to its transpose. The body-fixed thrust 


Equation (57) gives the thrust in the Z 0 direction. The 
first term (F £ 0 cos $) is the product of two small angles 
(£ and 0) times a sinusoid, which, when integrated into 
velocity, will average zero. Hence, the very good approxi- 
mation 

F' — F cos 0 (61) 

is made. 


Depending on how F' is used, cos 0 can or cannot 
equal unity. If the change in the velocity due to the effects 
of thrust misalignment and initial conditions is required, 
then the small angle approximation should not be em- 
ployed, However, if angular errors are of prime interest 
(as they are in this study), then little inaccuracy is 
incurred by letting cos 0 — 1. 


F. Velocity Vector 

Velocity can be determined from the value of the in- 
ertial thrust vector. The velocity vector is defined in the 
same way as the thrust vector, i.e., one component along 
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the inertial Z 0 axis, the other component defined in the 
complex cross plane. 


= 1, 

f K v dt 

mj 

f rV 

0 

rr 

/•' 

fc> 

/ Fldt 

m 

1 0 ‘ 


(62) 


then reducing, 


•-= 4 - }*■■«[•- >] 

-f r — ir g .Mi.x)i _ ii 

L (A+i) S „JL e l ] 


where m is the average mass of the body during thrusting , itN 

and g = 32.2 ft/sec-. If mass cannot be assumed constant, (1 4- A)$- 

then the integration of Eq. (62) must be done numerically. 

Equation (63) is the desired form of the cross velocity. 
Then, from Eqs. (60) and (62), Breaking into its coordinates reveals its character. 




— dt 


o' = = C i cos s^t + Cocoss u (l+A )t — (C,4-C 2 ) 

o' = 3m — C, sin sj 4- Co sin s 0 (1 4- A) t 4- C 3 1 


which is, when a t (Eq. 54) is substituted and integrated, 


t)’„ = - — (c ,4 «' - 1) 

m | s„ ’ 


rv 


iA 

s 0 (l f A) 


(e 




"I) 


where the C’s are the constants in Eq. (63). The above 
are seen to be modulated sinusoids, the first with a con- 
stant bias, — (Ci + Co), and the second with a ramp 
function (C 3 1) impressed over it. The general form of o' 
and o' is shown in Fig. 17. In the plots, either A ** —1 
and C 2 » Ci or A 4-1 and C 3 » C 2 . If A 0, the quali- 
tative solution shown is not valid. 


4- 


*0 ) 


Plots of o' vs c' are shown in Appendix B. These cor- 
roborate the cross-plotting of Fig. 17. 


where, temporarily, 


If cos 6 l, the integration of v’ z is straightforward. 
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However, if the small angle approximation cannot be 
made, i.e., the change in o' is desired, then 

o' = — f F cos Bdt = f cos \aAdt , (65) 

* ™Jo 1 1 

and numerical integration : s necessary. 


G. Velocity Dispersion 

It vjy = 0 due to the absence of initial conditions and 
thrust misalignments, then the resultant velocity would 
lie entirely upon the Z 0 axis. When 0, the tip of the 
resultant velocity vector is moved from the Z 0 axis (Fig. 
18) at an angle 



a = 



( 66 ) 


a is seen to be the velocity dispersion angle caused by 
the presence of a cross velocity, v ^ . The final value of a 
is reached when the time t — h, the burning time of the 
rocket motor. 


Equation (66) is expanded to (with the use of Eqs. 63 
and 64) 


«<«= t(v - 1 ) 


■t( 


N 


A(A+1) 2 s* (k+l)s, 


l)*o JV ) 


Fig. 18. Velocity dispersion angles 


and is treated as a random variable because of the ran- 
domness of the 0 O direction with respect to the true iner- 
tial system defined in Section VI. If it is required, it is 
seen from Fig, 18 to be tan ~'(v' y /v' x ). 

The wobble at the end of bum can be considered an 
output error. If this is denoted by then from Eqs. (43) 
and (44), 


1^1" (1+A)s„ 6 ° A{A + 1 )$l (68) 

T * is seen that if — 1 < A < 0, i.e., l z > 1 : , then the right 
term is negative and j fii I > 0 o . N itself is always defined 
positive (the worst case). 


+ 


iN 

U+aK 


(67) 


It is seen from Eq. (67) that if N = A (A + 1) s 0 0 o , 
the precession amplitude vanishes. The spinning ampli- 
tude can furthermore be nulled if the thrust misalignment 
{ = Oa (1 + X)/s 0 , but the bias term, N/( 1 + A) s 2 a , re- 
mains. Only if { = N = $ 0 — 0 does a(t) = 0 for all values 
oft. 

Of prime interest is the magnitude of a, i.e., | a | = 
(3(e i a + The radial direction of a with respect 

to the inertial system (angle r) in Fig. 18) is of no interest 


Another error which might be of interest (neglected 
here) is the error in the magnitude of the velocity vector 
(a is the error in direction). The value of o' is obtained 
frem Eq. (65); the odd term appearing in Eq. (57) (i.e., 
F „ 6 cos ^) still being neglected. This o' is compared to 
the value of the velocity had the thrust always been 
applied along the Z 0 axis, AOj = c in {m x /m f ), c being the 
exhaust velocity of the gases, m x the initial mass before 
thrusting and m f the final mass when t = ty, i.e., m; — m r 
= mass of propellant expelled. The difference in magni- 
tude of the two velocity vectors has been less than 1% 
for the range of parameters the author has dealt with. 
This means that the initial conditions and thrust misalign- 
ment cause less than 1% of the linear momentum of the 
rocket motor to be changed by spinning. 
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H. Limiting Cases 

A close inspection of Eq. (67) yields the following in- 
formation about a: 



recalling that N “ FiU I x . Hence, after t gets very 
large, a approaches a constant value which is in- 
versely proportional to spin-rate squared, and directly 
proportional to acceleration (F/m). This can be a 
useful approximation if the body goes through 
enough revolutions, so that the amplitudes of the 
first two terms in Eq. (67) are small compared to 
ctoc. If the parameters are further broken down, it 
can be seen that F = I T /t h where I T is the total im- 
pulse of the rocket motor, and h is the burning time. 
I T is a measure of the velocity requirement (since 
is a motor parameter). When the above value for 
thrust is substituted into Eq. (69), the result is 
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Figure 19 is a plot of a* vs $■ U for several values of 



. At the other extreme (as l grows small), after 
L’Hospital's rule is applied 

lim|a)= do = I £—0o I (71) 

t-> o 1 ' 

which is evident from the position of the thrust 
vector at t = 0 with respect to the inertial Z 0 axis. 

I. Taking limits on X 

As x — * 0 (sphere), a blows up in both the first 
two coefficients (Eq. 67), as would be expected, 
since a sphere has no body-fixed gyrodynamic sta- 
bility. However if t — > oc at the same time, the re- 
sults of Eq. (69) are 

, , . N 

inn f hm a | = — — ; 

x-o 1 / sj 

the same result is obtained if the limits are taken 
in reverse order and L’Hospital’s rule is applied. If 
t 0 as A 0, from Eq. (71), 

lim / lim a X = I £— 0 O I 
x-»o * o > 1 1 

As X-» —1 (thin rod), the second ana ~iird terms 
of Eq. (67) blow up and cause divergence of a. This 
also can be anticipated, since when X = — 1, l s = 0, 
and there is no angular momentum generated along 
the z-axis to counteract the effect of thrust mis- 
alignment. 

When X-» +1 (flat disc), the three coefficients in 
Eq, (67) are minimized with respect to X, thus min- 
imizing a ( t ). This results from the maximum mo- 
mentum developed in the z-dir^ction. 

1. Taking limits on s 0 

If s 0 0, a divergent form of Eq. (67) appears. 
This, however, is not indicative of the case since 
the small angle approximations no longer hold. Ap- 
pendix D treats this special case. 

As s 0 -» oo , the resultant error a($) -> 0, since in- 
finite stability is being approached. 


10 
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I. Maximum Envelope of Solutions 

If the constant coefficients in Eq. (67) are denoted K u 
K-, and o* respectively, i.e.,* 

a = 4 L ( C ” ,< “ 1) + — - 1) + ia* 

r v 

(72) 

the real and Imaginary components of a may be written 
as 

K 

a R — Qe a = cos s 0 t 

+ Aco SS .,1 + a ;1 -(M) 

K 

a / — <iJ/n ex = — ~ — sin Sq t 

+ -j- sin s 0 (l + A)t 4* a* (73) 

If the maximum values of the sinusoids are chosen in 
such a way that if both K, > 0, K 2 > 0, or both Kj < 0 
and K z < 0, (a w is always > 0), then 



Fig. 20. Maximum envelope of solutions 


The dashed lines are the real solution of | a(t) | in Eq. (67). 
It is to be noted that the body should go through at 
least 3 revolutions, i.e., s 0 ti/'2ir > 3 for this envelope 
approximation to be valid. 


a.<-=p(Ki + K*) 

«/<“-(K 1 + K 2 ) + o„ (74) 

or if Ki < 0 while K 2 > 0, then, 

*,<y(KI + kT 

« R < y(K 2 -K,) + a M (74a) 

The maximum magnitude of a will then be given by 

| a(t) | <(«•« + a*)* (75) 

If Eq. (74) is substituted into Eq. (75), it is seen that 
j a (t) | is a fourth order polynomial which approaches a 
hyperbola as a» -» 0. The curve at the point t = 0 has a 
singularity, and thus must be mated to the | a (f) j axis 
at Oo- Figure 20 shows how the envelope of | a(f) | behaves. 


tThe K’s used here are totally independent of those used in Sec- 
tion III, and are the result of an unfortunate choice of symbols. 


J. Some Examples 

A few examples are given of the results in this Sec- 
tion. The basic parameters will be the same as those of 
the last Section, i.e., Apollo and Ranger capsules, same 
spin-rate, etc. The error output of Section III is used as 
an input here (wobble angle 0 o ). Table 1 lists the perti- 
ment data. The results are shown in Figs. 21-27. 


Table 1. 

Basic parameters 



Apollo 

Ranger 

K * 

+ 216.4 

- 5,000, 
-10,000, 
-20,000 

Thrust misalignment f, rad 

0.004f 

0.004f 

Rocket burn time sec 

5.0 

9.6 

Capsule shape X 

Vs 


Initial spin rate, so, rad/sec 

1 2 5 10 

30 

Initial wobble angle B n , rad 

0.025 0.013 0.005 0.003 ** 

o.oor* 

*K — • « cotch-oll term arising out of Eq. | 67 ). 

A It 


** based on outputs of Section lit. 
tSoo Appendix C. 
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TIME, sec 

Pig. 22. a (f) vs t for Apollo-shaped capsule 
s 0 = 2 rad/sec, 0 O = 0.013 rad 



TIME, sec 


Fig. 24. a (f) vs t for Apollo-shaped capsule 
So = 10 rad/sec, 0 O — 0.003 rad 
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Fig. 25. a (t) vs t for Ranger capsule K — — 5000 


Three values of K are used for Ranger because of the 
large mass change involved during burning. In fact, the 
approximation for a ( t ) does not hold too well in ca^es 
such as this. 

Superimposed on the results of Eq. (67 ) are the maxi- 
mum envelope curves of Section IV, Part I. For the cose 
of s 0 = 1.0 rad/sec, d n = 0.025 rad (APOLLO), the enve- 
lope lies off the graph. As seen from the curves, the 
greater Lhe value of s 0 > the better the approximation. 



TIME, sec 

Fig. 26. a (f) vs f for Ranger capsule K = — 10,000 



V. STABILITY OF NOT-SO-RIGID SPINNING BODIES 


From classical rigid body dynamics it is known that a 
perfectly rigid rotating body is stable if it is spinning 
about the axis of either greatest or least inertia in torque- 
free space. If nonrigidity is present in the form of bend- 
ing, sloshing, rubbing, etc., then spin about any axis 
except that of greatest inertia is unstable. This insta- 
bility is in the form of a precession which causes fluctu- 
ating forces to dissipate energy in the nonrigid parts of 
-~3 body until the spin is about the axis of greatest inertia. 


When the Explorer I satellite was launched, the above 
phenomena was not taken into account, and the satellite 
(a cigar shaped body with long, flexible antennae ex- 
tending from the center, perpendicular to the axis of 
symmetry) was found to change its axis of spin from 
that of symmetry (least inertia) to almost a transverse 
axis (greatest inertia) in IV 2 hr. Subsequent analysis 
(Ref. 9) determined that the whip antennae served as 
ideal energy dissipators during precessional modes. 


.2 



JPL TECHNICAL, REPORT NO. 32-644 


Other possible applications of spinning for inertial 
stability have brought this question to light. As will be 
seen later, this p.nhlem is an extremely difficult one to 
solve if obvious sources of nonrigidity are not present. 
As an example, the whip antennae were removed from 
the satellite for the Explorer III launch, and it took 10 
days to v/obble 75 deg. This rate is slow enough that a 
number of items might have caused or contributed to it. 

This Section investigates the rate at which the angular 
velocity vector traverses from that of spin about an axis 
of minimum inertia to that of maximum inertia. The 
fluctuating acceleration on an element of mass is found, 
and a qualitative look at the effect of initial spin rate 
on wobble is made. 


A. Dynamical Equations 

Since the rotational energy is the coupling quantity 
between wobble angle 0 and time, 



where T is the rotational kinetic energy. dO/dT is easily 
derived from rigid body dynamics, whereas T is a more 
difficult term to find. The angular momentum vector re- 
mains invariant at J 0 regardless of the rotational energy 
variation (see Ref. 10). Figure 28 shows the position of J 0 
with respect to the body-fixed axes x, y, z, and 


where i, j and k are the body-fixed unit vectors. For a 
body of revolution about the z axis, 

li = 1/Wf , Jy 1/ Wy , J z — I z <J) Z 

and the rotational energy is 

Jl 

(77) 




!)+■ 


2L 21, 


and from Fig. 28: 


J5 + /* = /; sin' 6 

/*=/; cos** (78) 

where J 0 = l : s„ and is constant. Equations (77) upon 
substitution and rearranging, become 


2T _ sin 2 0 cos 2 9 _ 1 /I 

Jl h h L \h 



cos 2 9 


If energy is dissipated, 6 must increase as T decreases, 
which implies that l x > l~ for a real solution. This con- 
firms what was said before, that the system is tending 
toward the axis of highest inertia. 


If noth sides of Eqs. (77) are divided by the initial 
rotational energy (T 0 = 1^/2), and the ratio T/T 0 is 
defined as T r , then 


T r = 1 + X sin 2 9 (79) 


Jo = JA + Jvi + Jzk 



Fig. 28. Resolution of angular momentum vector 


where X = I*/l x — 1. This equation is of marked interest 
in that it states that the way in which the kinetic energy 
decreases is only dependent upon the shape of the body 
(X) and the wobble angle (9). This fact will be used later 
in determining the effect of initial spin on 9. Equation 
(79) is plotted on Fig. 29. 

The rate at which 9 changes with respect to T r is 

d9 rrfTVT* _ 1 

dT r " j_ d9 J X sin 29 

1 
2 

which is shown in Fig. 30. It is seen that d9/dT r has no 
real solution outside the limits (1 + X) < T r < 1. These 
bounds, of course, define the limits on the wobble angle 9, 
i.e., T r — 1 + X is total stability 9 = 90 deg, and T, = 1 
is the condition when 9 = 0. 


{(Tr- 1) (X+ I - T r )} 


(80) 
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x 

HO -0 9 -C 8 -07 -06 -0.5 -0.4 -0.3 -0 2 -0.1 0 



ENERGY RATIO, Tr 

Fig. 29. Kinetic energy ratio vs wobble ongle 



ENERGY RATIO, T r 

Fig. 30. Range of solutions for 66/61 R 


It is important to note that both the function 6 = f(T r ) 
and its derivative 6' *= d6/dT r are asymptotic at the 
limits of T f . Therefore, neither condition is ever really 
reached. If 6 = 0, i.e., T r — 1, all acceleration is normal 
to the spin axis and no fluctuating forces are present. 
Hence, the body is stable. By the same token, as 6 -» 90 
deg, the accelerations are approaching a constant magnb 
tude from the other direction, which indicates that 6 = 90 
deg is never reached, only approached. Therefore, when 
the wobble angle is said io be going from 0-90 deg, what 


is really meant is that it goes from some initial 12 value 0„ 
to some other value, which, in the case of Explorer ///, 
was 75 deg for the 10-day number quoted. 


The value of dd/dT , which is needed in Eq. (76), is 
then 


dLV r dT r y 1 

dT L 0 dO T 0 A sin 26 


(81) 


which is seen to be an inverse square function of the 
initial spin s 2 . 


B. Internal Forces 

The next question of interest is with regard to the 
forces which cause the internal bending, sloshing, etc. 
Initially, each element of mass dm in the body has im- 
pressed upon it a large centrifugal acceleration i p 2 p (see 
Fig. 31) which gradually damps out. The fluctuations 
arise out of the coupling of ^ with <£, the precession rate. 
For this portion of the analysis a new coordinate system 
is defined, also based on conventional Euler angles. ' 

Let an orthogonal triple £, 77, z be defined along the 
body-fixed roll axis, the origin being z above the inertial 


12 Or as an input, 0, from Section IV, which is the final wobble angle 
at the end of the thrusting phase. 



Fig. 31. Rotating coordinate system 
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origin The £ axis is constrained to pass through the Z u 
axis at all times. Hence, the new axis rotates about the 
inertial X.rY.rZ,, axes at a rate <£. The distance 2 defines 
the distance along z where the element of mass dm is 
situated. The fact that 0^0 will be neglected here, 
since the rate of change of wobble angle is normally 
orders of magnitude less than 6 and p. It is also seen 
that the discussion of the forces on dm will be inde- 
pendent of $ , hence 6 may be explicitly omitted, since 
the coordinate system will be defined when <j> — 0. 

The problem may then be stated as: given a coordinate 
system (£, y, z) rotating about an inertial system 
at an angular rate <j> which will be time variable. A par- 
ticle dm at a distance p from the .v axis is rotating in the 
£, 7/, z system at a rate \p. What are the nonconstant forces 
on dm? 

Since p and z are constant with time, in the general 
equation for the acceleration of a point in a rotating 
reference frame (see, for example. Ref. 11, p. 210), all 
but the centrifugal and coriolis terms vanish. The mass 
dm is seen to have a centrifugal acceleration p ^ with 
respect to the moving z-axis; and as seen in Fig. 31, an 
acceleration <j> 2 (z sin 6 — p cos 0 cos \p) with respect to 
the Z 0 -axis. Since the fluctuating terms are of principal 
concern, one pertinent force parallel to the X, -Y„ plane 
plane is 

F j = -<#>- p cos 8 cos p dm (82) 

The coriolis term in the general acceleration equation 
is 2 u) X v, where co is the $ vector as reflected in the 
£, rp z axes. 

9 9 

a) = <f> sin 8 i 4 </> cos 8 k 


F , and F-, can be vectorially added to determine that the 
component of force normal to the £-7/ plane (i.e., along 
the --axis) is 

F r ~ (2^ <f> - <f> 2 cos 8) p sin 8 cos $ dm (84) 
and the component in the £-7/ plane is 

F t . v — F } cos 8 -- — <fy“ p cos vif cos- 8 dm 

(85) 

However, the important thing to note is that the force 
in either plane can be expressed as 

F = / (p, <£, 6) cos $ dm 

C. Variation of 4> and 1 ft with 0 

As 6 increases due to the internal dissipation of energy 
and the body seeks a higher inertia, spin is transferred 
from ^ to ^ in such a way as to conserve the angular 
momentum. The relation between these two is found from 

J z ~ L ~ ] t) cos 8 

and - Z : So where s„ is the design spin 

rate, ideally all about the z-axis, and 

• • 

CO; “ COS 6 ~f- \j/ 

From the first two equations above, = f (8) is deter- 
mined where substitution in the third yields 

ip ~ (s„ - j>) cos 8 (86) 

• • 

Thus, it is seen that when 8 = 0, \p 4- $ — The de- 
pendence of 4> upon 0 is seen from 


where the unit triple i, }, k is defined in the £, yp z frame. 
The velocity of dm in the £. 77, z frame is 

v p \p sin $ i -t- p ip cos \p j 

and the coriolis acceleration is then 

* * 9 * 

2 (o) X r) = 2( — p\p<t> cos 6 cos $ i-l p $ <f> cos 8 sin \p j 
+ p \p </> sin 8 cos \jj k 

The magnitude of this in the £-77 plane is seen to be 
constant, p\p<p cos 8 rotating at a frequency Thus, an 
observer sitting in dm would not experience a fluctuating 
force in this plane. The fluctuating force is along the 
%-axis and is 

F z = 2 p \p <j> sin 6 cos $ dm (83) 


n + Ji = iite+»D = n™** 

-= <j> sin 8 sin ip 0 cos \p 
0 ~ 4> sin 6 cos ^ — 0 sin \p 


whence it may be determined that 

<£ 2 ” s 2 (A -f 1)' J — 6 - esc 2 6 (87) 

Here it is seen that if 0 is assumed zero, then <f> is con- 
stant and the device for transferring the spin from p and 
<j> is no longer real. When 0 < 8 < 90 deg, <£ is well 
defined, and as 8 90 deg, the right hand term becomes 

8~ esc 8 &- and 


4>' 2 - + ~ ( 88 ) 
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for a final value. As ft — » 0, the <j> and ^ tend toward coin- 
cidence, and the orthogonal system loses one degree of 
freedom and the solution of Eq. (50) has as a root sin ft 0. 
Since at ft — 0, w v — <J> i- the decree of freedom may 
be used by stating that — 0. This assumption ’s 
apparent since in Eq. (82), F, must vanish when ft = 0 
(all acceleration is constant), thus <£ = 0. In the same 
manner, F> vanishes at ft ~ 0: when ft — 90 deg, all accel- 
eration is along the z- axis, hence F^ ~ 0; and since ^ = 0 
(Eq. 86), F 2 = 0. 

Defining = 0 when ft = 0 and <j> = s 2 0 (A + l) 2 — 

ft 2 esc 2 ft , the moment <j> > 0 might seem inconsistent. 

However, the coupling between <£ and \p (through ft) is 
perfectly valid when ft > 0. In fact, Eq. (87) is exact. 
It is at the singularity ft — 0 that difficulties arise. Since 
<£ and \p are indistinguishable then (with respect to the 
mathematical interpretation, but not to the physical), 
4> may arbitrarily be picked. 

The forces on the element dm are then seen from Eqs. 
(84) and (85) to be (assuming 0 — 0 in Eq. (87) 

F, = — $ 2 0 (1 4- A) p sin ft cos ft cos \p dm 

F cv = — s 2 0 (1 + A) 2 p cos 2 ft cos \ff dm (89) 

• * 

The value of ft to be used in determining <£ (Eq. 87) is 

that in Eq. (76), which in turn needs dT/dt as an input; 
hence, an iterative solution is required. If, at the end of 
spin-up, injection, etc., and ft l} are known, then <£ n can 
be determined from Eq. (51). Ft and F 2 for the first in- 
crement of time can be found. Then from experimental 
damping data, energy dissipation is determined, i.e., 
dT/dt. ft is calculated from Eq. (76) and substituted into 
Eq. (87), where a new, larger value of is found. The 
iterative loop is continued as ft 90 deg. 


Let the motion of the element of mass dm be approxi- 
mated by a linear oscillator, the equation of mobon of 
which is 14 

dm^li + it + o>l = F 1( cos i t (90) 

where n>„ is the natural frequency of dm , a the displace- 
ment, and the ratio of total energy to energy dissipated 
in one period. F„ is the magnitude of the fluctuating driv- 
ing force discussed in Section B, and can be characterized 
as F n = s 2 C, (A, ft ), where C, is independent of s 0 , and 
is only a function of the shape of the body and the 
wobble angle ft. The frequency ij/ can be written as 
s C 2 (A, ft) from Eqs. (86) and (8 7 ) (letting ft = 0). The 
solution can then be written as 

5 2 F cos (it 4-8) 

u = .. , ( 91 ) 

VK - cj »; )* + (».. cttjqt 

where 

Wo C > Sn 

8 = tan “ “ ()(«; -c* Jf) 

Note here that the time varying value of F n was inte- 
grated as though it were a constant. This is allowable 
for two reasons: first, $.> is usually orders of magnitude 
greater than ft , the driving force behind F n ; and second, 
an attempt is being made to determine the effect of vary- 
ing s n on identical bodies (i.e., A’s are equal). This 
implies that the way in which the body varies its internal 
energy is determined by Eq, (79), which is a function of 
ft alone once A is stipulated. Hence, for comparison pur- 
poses, it little matters what form F„ assumes. 

The energy dissipation is obtained from the energy 
integral when nonconservative forces are involved and is 

(92) 


D. Effect of Initial Spin on Tumble Rate 13 

Since the ratio of the kinetic energy T to T 0 is inde- 
pendent of initial spin (s 0 ) the motion of the same body 
(A) can be compared at several different s 0 . This is be- 
cause the motion will always be of the form in Eq. (79), 
dependent only on A and ft. 


13 This Section is an expansion of a private communication from 
Professor Leverett Davis, Jr., of California Institute of Tech- 
nology to the author. 


where Ff } is the nonconservative force, which, in this case, 
is the damping term, 

£?, - ^itdm (93) 

whence 

dT wo • 2 i 
1t-Q uim 


14 The equations and analysis here can for the most part be found in 
the first two chapters of Ref. 10. 
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u is calculated from Eq. (91). Since u" is an oscillating 
function [sin- (ipt + 8)J, the average value of T over 
one cycle is 


/dl 

\dt 




dm 


s G C 2 F 2 
^0 U 2 r 0 


(-j -c*. »:)■+(— y ! 


0J 0 J 

2 Q dm 


(94) 


The equation for 0 is then 



— “ 


r i i 

de _ 

C 2 F 2 
^2 0 


Q 0,0 S 0 

dt 

2 h A sin 2 6 


[« ^ * 2 . y + ( u,,,c q ") J 


(95) 


Three case^ present themselves: 


e. When the natural frequency of the system is low, 
i.e., s 0 > > w 0 j thm 


t 

and 0 (and time) is independent of $ 0 . 

If p passes through a resonant frequency, = w 0 
at some point of the transition, then the energy 
dissipation will be greater. The magnitude of 0 is 
also a strong function of w 0 and Q. In case a, where 
w 0 > > s 0 > > s 0 /Q , then 0 a < P. In case b, where Q 
is small, and w 0 and s 0 are intermediate, it may be 
approximated that 0 b * - P < 1). In case c, where 
w 0 is small and O is in f ermediate, Cl < < 1, and it 
follows that 9 C > P. It may be generalized that 

e c > o a 

Explorer I, with its long whip antennae with low 
natural frequency, fell in case c. 

As a rough cut of the time required to tumble, if 
's^dT/dtyis assumed constant over the range of interest, 
then 


a. If the natural frequency is very high s gid system), 
so that w 0 > > s 0 , and the damping is very small, 
Q >> 1, which implies w 0 >> $ 0 /Q , then Eq. (95) 
reduces to 


M = 


Tmin 


<dr/d#>a 



where ^dT/dt^^y is determined from Eq. (94). 


where 

C 2 F 2 

p = . . r jL, 

2 I s A sin 26 

and the rate of change of 6 is proportional to 
s 4 0 . Hence, the time required to go the same number 
of degrees (say from 2 to 75) would be proportional 
to i/«; . 

b. If Q is very small (large damping), so that 
So/Q > > to 0 > > then Eq. (95) becomes 


£. Discussion 

The problem may be reduced to this: Since the dynami- 
cal relations (Section V-A) and forces on the element 
(Section V-B) are known, it remains only to (1) recognize 
those elements of the body which will be affected by the 
forces, (2) define the mathematical model to which each 
belongs, and (3) determine the values of the constants to 
be used. 

General conclusions are best summed up by describing 
items to be avoided when designing bodies for spin. 

1. Try to make roll inertia greater than the inertia in 
other axes (e.g., Apollo), in which case the problem 
vanishes, 


B h = P 


Q*l 


Wo 


Cl 


2. If this cannot be accomplished, make the body as 
rigid as possible, avoiding any liquids, dangling wires, 
etc. Make all bulkheads parallel to the roll axis stiff. 


and the time to go between equal values of 6 is 
proportional to 1 /s 2 0 . 


3. If neither 1 nor 2 is possible, keep the design spin 
rate s 0 as small as possible within other system constraints. 
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In general, it appears from past results that this phe- 
nomenon is not hard to design around. If, by merely 
removing the whip antennae on an Explorer Z, the tumble 


rate was decreased by over two orders of magnitude, then 
one may conclude that internal damping is something one 
designs in rather than out of the typical spinning body. 


VI. SUMMARY AND GENERAL CONCLUSIONS 


A . Error Plan Summary 

The overall errors discussed in the Report are depicted 
in Fig. 32, and can be somewhat broken down into chro- 
nological order. 

Section III 

a. Initial errors of spin-up. These are two: a tip-off 
angular velocity 0 O , and an error angle at time of 
spin-up which is 0 O + 8 0 St , 0„ being the initial dis- 
placement, being the time from release to initia- 
tion of spin-up. 


b. Dispersion of the angular momentum vector from 
the inertial axes during spin-up (fl/). 

c. The wobble angle produced during spin-up (/?). 

Section IV 

a. Velocity vector dispersion angle (a) produced by the 
initial conditions and the misalignment of the thrust 
vector. 



TIME -► 

Fig. 32. Overall error plan 
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b. Velocity vector dispersion magnitude. 

c. Final wobble angle (/?,). 

Section V 

a. Effects of nonrigidity on wobble angle. This part is 
not an error (as in Sections III and IV), however 
it can affect the spacecraft performance if neglected. 


B. Genera/ Conc/us/on$ 

The first two conclusions may be stated categorically, 
(1) whenever possible, A > 0, i.e., the inertia about the 
spin axis should be as great as possible. If this is true, 


errors in Sections III and IV are minimized with respect 
to shape, and the discussion in Section V does not apply; 
and (2) all error sources, such as initial conditions, thrust 
misalignment, etc., should be minimized. Other conclu- 
sions are: (3) during spin-up, both the spin torque and 
the spin rate should be as large as possible; (4) when 
thrusting, the number of revolutions through which the 
body turns should be as high as possible, i.e., s tJ h should 
be large. As a rule of thumb (noted from the results of 
Section IV), if sj b /2w < 3 revolutions, then little is gained 
by spinning; and (5) for coast stability, the body should be 
as rigid as possible and be spinning at the lowest possible 
spin rate. The obvious contradiction between conclusions 
4 and 5 is resolved by despinning after thrusting. Appen- 
dix E contains general equations for despinning by the 
4 yo-yo” method. 
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APPELulX A 

Orthogonal Transformation for Analog Simulation 


If an analog simulation of spinning dynamics is desired 
(including the angular rate transformation from body- 
fixed to inertial coordinates), the conventional Euler 
angles defined in Section II are found to be unsatisfactory. 
This can be seen from Eq. (50), which gives the angular 
rate transformation 


<a T 


“sin 6 sin i p 

COS \ff 

O' 

0) v 

= 

sin 0 cos \f/ 

^sin \ff 

0 

01* 


cos 0 

0 

1 


For the solution, the inverse of the above equation is 
required, or, if the matrix is A, 


1 

i 

< 

< 

II 

1 

e & 
1 

L*J 

’ ‘ 1 

L**J 


where | A | is the determinant of A and j | A j j is the adjoint 
of A. By inspection, j A | = — sin 9, and since 9 is nomi- 
nally zero, or if not, almost zero, the equations in (A-l) 
are unstable for analog simulation. Another set of Euler 
angles must be defined in which the determinant of the 
angular rate matrix is not nominally zero. Figure A-l 
shows one possible transformation which is stable over 
the region of interest. The rotations, in order, are 



and the angular rate transformation is 


[A] = 


"1 

0 

0 


0 0 " 
cos A sin A 

— sin A cos A 


[G] = 


cos G 
— sin G 
0 


sinG 
cos G 
0 


O' 

0 

1 


[ 0 ] = 


|"cos0 

0 


sin 9 


0 — sin 9 * 

1 0 

0 cos 9_ 


<o x 

<Oy 

<»* 


cos G cos 9 — sin 0 

— sin G 0 
cos C sin 9 cos 0 



(A-2) 


which can be verified from Fig. A-l. The determinant of 
the matrix in Eq. (A-2) is — cos G, which if the body-fixed 
y-axis is defined as that of spin, is nominally unity. Hence, 
stable solutions of the inverse of Eq. (A-2) result. The 
fully expanded transformation of the three rotations is 


[0] [G] [A] = 


cos G cos 9 

— sin G cos 9 + cos G sin A sin 9 
cos A sin 9 


sin G cos A 
cos G cos A 
—sin A 


— cos G sin 9 + sin G sin A cos 9~ 
sin G sin 6 + cos G sin A cos 9 
cos A cos 9 


whence any body-fixed vector v is transformed into an 
inertial vector F 0 by 

Vo= [A]' 1 • [G]- 1 • [0]- 1 v 
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APPENDIX B 

Exact Formulation and Numerical Solution 
I. NONLINEAR EFFECTS 


If nonlinearities are present in the system, die basic 
differential Eq. (41) and subsequent equations cannot be 
explicitly solved. These nonlinearities usually arise out of 
more exact formulations of the problem. 

a > Variable Mass During Burning. The largest single 
cause of inaccuracies in Section IV is the mass decrease 
during burning, which decreases inertia. A good mass 
approximation assumes a constant thrust, thus constant 
mass flow rate, 


m 


= m 0 - mt = m 0 - (Jr^j 1 ( B ‘l) 


where m 0 is the total mass at t = 0, and I S p is the rocket 
propellant specific impulse. Inertia changes may then be 
characterized. If the rocket motor is considered a point 
mass at the end of a thrusted body (see Fig. B-l), then 
the inertia change in the cross axes (pitch-yaw) is the 
initial inertia decreased by a changing mass effect as 
determined by the parallel axes theorem. In effect, the 
“rigid body” mass ( M ) remains constant, while the rocket 
motor attached to one end decreases in mass (m). The 
pitch inertia (as a function of time) is 

h(t) = I u + Im 4- fifx 2 + m (l - x) 2 (B-2) 

where Iu = inertia of M (rigid body) about its c.m. 

I m = inertia of m (rocket motor) about its c.m. 

I = distance between mass centers (can be 
variable). 
m 


x = 


M + m 


l 


ROCKET MOTOR c.m. 

-H p h 


POINT OF THRUST 
APPLICATION 



-COMBINED c.m. 

M- *-H 



RIGID BODY c.m. - 

Fig. B-l. Schematic of mast variation 


and m(t) is given by Eq. (B-l). The torque about the 
pitch axis is then 

L, = Ft(p + l~x) (B-3) 

If the assumption of m made in Eq. (B-l) is not valid, then 
whatever m(£) is assumed must be pointwise integrated. 

Decreases in the roll inertia can be neglected, since to 
a first order the mass thrown out by the motor takes with 
it angular momentum, keeping the spin rate about the 
roll axis constant at Except in cases of very large 
mass decreases during spinning, the variation of I- can 
probably be neglected (since the mass of propellant being 
expelled is usually situated close to the roll axis, thus 
contributing only a small part of the total I z . In the case 
for Apollo, all changes in inertia, both pitch and roll, are 
neglected. However, for Ranger, when % of its mass is 
expended during burning, these variations must be taken 
into account. 

b. Jet Damping . When a body, rigid or not, has a 

cross-angular velocity (i.e., = o> x -f iw v ^ 0) during the 

expulsion of mass along the z-axis, a damping arises out 
of the fact that the cross causes a “wagging” of the jet 
stream, i.e., a change in the direction of the linear momen- 
tum vector. The magnitude of this damping (Ref. 8, p. 21) 
is proportional to the rate of mass expulsion, and to the 
cross-angular rate. The damping occurs because, if <o 
increases, the rate of change of the linear momentum 
vector of the exhaust gases is greater, and thus causes a 
greater resistance to the increase in ai. If m is the mass 
flow rate, and the distance (l + p — x) from the combined 
c.m. to the point of thrust application, the jet damping 
term is — w m (l + p — x) 2 . 

c. Small Angle Approximations . The expediency of 

letting sin 0^6 and cos 6 1 in the inertial transfor- 

mation is obvious from the analysis. From practical 
considerations, this approximation is usually valid, since 
real systems seldom can tolerate angular disturbances 
greater than 5-10 deg without intolerable error buildup. 
If, however, this approximation cannot be used, the 
solutions of the inertial transformation must be made 
with the full Eulerian expansion, resulting in transcen- 
dental differential equations. 
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d. Roll Moment During the burning of some rocket 
motors, a roll torque is generated. If this is to be con- 
sidered, the w r vector in Eq. (41) is no longer constant, 
but must be varied according to the roll acceleration. 

e. External Forces . Gravity and solar pressure are in- 
ertial forces, hence they may be neglected with respect 
to the spinning dynamics, but must be included in any 
velocity calculation (e.g., Eq. 62). Aerodynamic forces 
are only partially inertial since they depend on the angle 
of attack 8 (sometimes on <j>), and are independent of ^ 
in a body of revolution. Good discussions of aerodynamic 
forces can be found in Refs. 4 and 8. 

f . Cross Products of Inertia. It has been assumed in this 
report that all cross-products vanish, since the present 
technology is such that the principle axes of the body 
and the geometrical axes can be made almost coincident 
in bodies of revolution. To consider the cross products, 


the full inertia tensor must be included in Eulers dy- 
namical equations and carried along in the subsequent 
analysis (see, for example, Ref. 1 or 12). 

g. Nonrigid Effects. If the strain on structural mem- 
bers and mass deflections are to be considered, the 
problem inert.^es tremendously. Euler’s dynamical equa- 
tions must be modified to include cross-products of inertia, 
time-varying coefficients, and ether nonlinear effects 
discussed below. Each effect of nonrigidity must be 
treated separately, e.g., mass unbalance will affect inertia, 
cause cross-products, c.m. shift, etc.; straining of the 
structural members will cause mass unbalance, rotation 
of the thrust vector, etc. 

h. General Nonlinear Effects . These include variations 
in thrust vector misalignment with time, variation of 
thrust with time, etc. Each one of these must be included 
in th? basic differential equation and be numerically 
integrated. 


II GENERAL FORMULATION 


Combining as many of the aforementioned effects as 
possible, the basic equation for the body-fixed w is then 

i> + i u„(t) { A(f) - m(l + p-x)-)<u = y|j- (B-4) 

where the time functions of <•>-, A, x, l, L, /, and m (if not 
constant) are to be stipulated. Equation (B-4) is then 
numerically integrated to determine the cross-angular 


rates w(t) = <■>,(*) + i <u„(t), which, with the solution of 
, 0 - = N : (t) [if there is a torque in the ^-direction] are 
substituted into Eq. (8). The numerical solution of this 
set of equations will give the time functions of the Euler 
angles (<f>, 0, <fi) which define the inertial transformation 
matrix T. The inertial thrust vector is then determined, 
which, when divided by m(t) and integrated, yields the 
inertial velocity vector. 
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III. SPECIFIC EXAMPLES 


Ranger is an ideally suited example here, since % of 
the weight of the original system is expended. Table B-l 
gives the parameters used in the problem. 


Table B-1. Ranger parameters 


Thrust F (constant), lb 

5500 

Burning time f*, sec 

9.6 

Payload mass, lb 

96.0 

Flow rate M (constant), Ib/sec 

20.0 

Spin rate so, rad /sec 

30.0 

Lever arm in (!) (variable), in. 

See Fig. B-2 

Pitch inertia /„ lb-in. 2 

See Fig. B-2 

Inertia ratio, X 

See Fig. B-2 

Capsule mass, M, lb 

See Fig. B-2 

Thrust misalignment rad 

0.004 

Initial wobble angle rad 

0.001 

H FI 


K = TJ7 

See Fig. B-2 


This problem is identical to that for Ranger in Section 
IV-J, except there, three points ( K = —5000, —10,003, 
— 20,000) were used to examine the effect of K on a. 
These three points are seen to represent the K vs t curve 


on Fig. B-2. Looking at Fig. B-3, the output of the above 
problem, it is seen that the K -- — 10,000 in the approxi- 
mation is close to the true solution on Fig. B-3. However, 
the shape of the solutions do not agree too closely. If 
the proper value of K could be assumed, the approximate 
solution could suffice for values of a(t). 

The Apollo problems in Section IV-J were also run on 
the big program V> illustrate two points. If the mass 
change is small, * e a(t) in Figs. 22-25 do in fact closely 
approximate tf • ’.rue solutions in Figs. B-4 through B-7. 
The other rea^ s that the larger program allows a look 
at V' vs V' y 1 a f real vs a f imaginary. The Apollo 
problems were hose ", because the high spin rate of the 
Ranger capsule made the . oo confusing. Figures 

B-8 through B-li show V* vs V }/ for the problems listed 
in Section IV-J. The shape ot the curves corroborates the 
statements made in Section IV-F. Figures B-12 through 
B-15 show a, real vs a f imaginary. This is the trace that 
a ray of light located along the body-fixed z-axis would 
make as seen by an observer sitting above the capsule in 
the inertial frame. 






|a(r)J,rod „ la (/)!, rod 
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Fig. B-5. a U) vs t for Apollo-shaped capsule 
(vuriable mass), s 0 = 2 rad/sec, 0 o = 0.013 rad 



Fig. B-6. a (f) vs f for Apollo-shaped capsule (variable 
mass), s 0 = 5 rad/sec, 0 o = 0.005 rad 



Fig. B-7. a it) vs t for Apollo-shaped capsule (variable 
mass), s 0 = 10 rad/sec, 0 o — 0.00? /ad 



Fig. B-8. V* vs V y for Apollo-shaped capsule, 
So = 1 rad/sec, 6 0 = 0.025 rad 
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Fig. B-13. 3le o, vs Sm 0 ; for Apollo-shaped capsule, Fig. B-14. £l{e a, vs Sm a, for Apollo-shaped capsule, 
*o = 2 rad/tec, $ 0 = 0.013 rad s 0 = 5 rad/ sec, $ 0 = 0.005 rad 



Fig. B-15. a / vs Sm a t for Apollo-shaped capsule, 

So = 10 rad/sec, 8 a = 0.003 rad 
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APPENDIX C 

Thrust Misalignment of Solid Propellant Rockets 


'9 An attempt is made in this appendix to justify the use 
of £ = 0.004 rad for effective thrust misalignment from 
data gathered from the in-flight measurements of the 
second and third stages of the Scout launch vehicle. All 
four stages of the Scout are solid propellant rockets; 
however, the first stage burns in the atmosphere, making 
any misalignment data almost impossible to separate 
-zf from the aerodynamic effects, and the fourth stage is 
spin-stabilized. The middle two stages are attitude stabil- 
ized by bang-bang peroxide jets in pitch/yaw/ roll.* 
Vehicle pitch and yaw motion and rates are telemetered 
and recorded. The slope of the pitch rate trace when the 
H 2 0 2 jets are off gives the vehicle acceleration in the 
pitch plane; yaw acceleration is obtained in a like man- 
ner. The vehicle nominal thrust, inertia, mass, etc., is 
known as a function of flight time. Hence, the effective 
thrust misalignment may be determined in either the 
pitch or the yaw plane.** Figures C-l and C-2 illustrate 
how the thrust misalignment varied on Scout flight S-127. 
Table C-l lists m?Yimum thrust misalignment in degrees 
for the second and third stages for 12 flights, £ p being the 

♦Roll deadband is ~ ± l Vo d minimizing pitch /yaw coupling. 

♦♦This thrust r.dsalignmtnt is effective because it is what the vehicle 
sees from r torqumg standpoint, hence it includes the effect of 
center-of-mass variations, vehicle bending, etc. It is therefore 
erroneous to speak if as only a solid rocket thrust misalignment. 



Fig. C-2. NASA Scout S- 1 27 3rd stage pitch and yaw 
thrust misalignment vs time 

component in the pitch direction and being the com- 
ponent in the yaw direction. The maximum values do not 
necessarily occur simultaneously. 

As seen in the table, the value for £ rarely goes above 
0.2 deg for stage 2 and stage 3. The weight of peroxide 
loaded on board for these two stages is based on a 0.25- 
deg thrust misalignment for the full duration of second 
stage bum, and 0.10-deg thrust misalignment for full 
duration for third stage burn. The peroxide supply has 
never been exhausted. Thus, from the above numbers, 
it is thought that a value of £ of 0.004 rad (0.23 deg) is 
considered conservative for the examples considered hi 
this Report. 


* 



Fig. C-l. NASA Scout S-127 2nd stage pitch and yaw 
thrust misalignment vs time 


Table C-l. Scout thrust misalignment 


Vehicle 

Sotond itage 

Third stage 

number 

max, dog 

max, dog 

£ p max, deg 

U deg 

ST-6 

■ 

0.100 

0 074 

0.026 

0.011 

ST -9 

0.228 

0.048 

- 

- 

S-U3 

C.079 


t 

t 

$-114 

0.096 

0.064 

- 

- 

$-115 

o.i ;? 

0.136 

*.048 

0.077 

$-116 

0.090 

0.012 

t 

0.020 

$-118 

0,127 

0.02/ 

0.071 

0.182 

$-119 

0.215 

0.103 

0,065 

0.034 

S-120 

0.134 

0.015 

0.1 09t 

0.455ft 

$•122 

0.159 

0.029 

0.069 

0.092 

S-127 

0.130 

0.014 

0.062 

0.077 

S-132 

0.177 

0.077 

0.110 

0.016 

t&ody otd licit! on* obscured Hi* “mean" 
ft Low thrust prior to burnout. 

rate data. 
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APPENDIX D 

ZERO SPIN RATE 


In certain applications, a short, high thrust with no 
spin may produce acceptable velocity dispersions (if the 
velocity increment is low enough). If such is the case, 
the small angle approximation made throughout Section 
IV no longer holds. Additionally, if the thrust misalign- 
ment is assumed constant (as it has been), the motion 
when $ 0 = 0 is confined to two dimensions; either X 0 or 
Y 0 ".nd Z 0 . Referring to Fig. D-l, the curved path is the 
actual velocity path followed in the X 0 -Z 0 plane, and the 
straight line (V) is the resultant velocity vector. 

The angular acceleration of the body is a function of 
the torque, i.e.. 


8 - K = ^ 

l x lx 


(D-l) 


A) 



Fig. D-l. Velocity diagram when S 0 = 0 


whereupon, if the initial tip-off rate is 6 0y 
• * 

6 = 6 o H — ~ t, and integrating once again, 

lx 

0 = 6 0 + 0„ t + (D-2) 

in Section III, 0 C — 0 through the rotation of the 
inertial axes an angle 0 O . The velocity components are 
then seen to be 


The integrals in Eq. (D-3) are seen to be modifications 
of the Fresnel integral, and if = 0, 



V, 




sin $ dt 


/.' si 

= — f ( sin 0 O t cos -7^- 1- + cos d„ t sin t*\ dt 
m Jo \ 27, 21, ) 

1 V, = i£f co sOdt 

m Jo 

= [ ( cos 0„ t cos 1- — sin 0 O t sin tA dt 

w Jo \ 2I X J 

(D-3) 


where C and S are defined in Section III-F. The error in 
the velocity direction is then 


a = 



(D-5) 


and the error in the magnitude is ^1 — 

Unless tb is very small (< 1 second), V, will probably 
be excessive. The velocity is necessarily kept small so 
that the thrust misalignment is reasonably bounded. 



.4 
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APPENDIX E 

DESPINNING 


After a body has been spun-up and a maneuver has 
occurred, it is often necessary to de-spin. The following 
is a brief listing of some despin devices: 


for the rigid cord yo-yo are given below, and are taken 
from Ref. 14. The deployment of a rigid cord yo-yo is 
depicted in Fig. E-l. 


1. Gas jet (hot or cold). The main drawback in addi- 
tion to high weight is the fixed total impulse, which 
will despin only a fixed number of rpm. That is, if 
the jets work perfectly, a 10-rpm error in initial 
spin rate will also appear in the final spin rate. 

2. Body-fixed magnetic rods in a magnetic, field can 
also stop a rotating body (Ref. 13). However, this 
method has four main drawbacks: a magnetic field 
such as the Earth's) must be available, the time 
involved to stop the rotating body is considerable, 
all spin must be removed, and the use of magnetic 
rods in the body might have some interactions with 
other elements. 

3. “Yo-yo s” are devices which decrease spin by re- 
leasing weights on cords which unwind due to the 
spin, causing large increases in inertia about the 
spin-axis. When the cords have unwound, they are 
released, taking with them the difference between 
the initial and final values of angular momentum. 

Yo-yo s with both rigid and stretch cords have been 
investigated extensively in the past five years (Refs. 9, 
14, 15, 16, 17). For reference, the necessary equations 


If <o„ = initial spin rate, rad/sec 

6 = final spin rate, rad/sec 

y — Q/o> 0 

m ~ mass of despin weight (for 1 cord), slug 

I ~ inertia of spinning body about spin axis, 
slug-fit 3 

a -radius of body (to cord attachment), ft 
R =- required length of cord, ft 
T, 1U , = maximum tension in cord, lb 


t = time to despin, sec 


then the length of one cord (there are 2 as shown in 
Fig. E-l) is 


R = 


ni-y) (r-r2ma-r p 

L 2m(l-ry) 


(E-l) 


Note that this equation states that for a given R and m, 
y = 0/ a>« is constant. Hence, yo-yo’s take out a given 



(a) 



Fig. E-1. Rigid cord “yo 
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proportion of spin, and a 1% error in initial spin iate 
will remain a 1% error after despin. However, if 6 — 0, 
then Eq. (E-l) above becomes 



which states that to despin to 0 rpm, the dimensions of 
the despin mechanism (i.e., R. t and m) are independent 
of Wg. The time to despin is found bv letting ft 
in Eq. (E-l) and solving. 


, 1 f (1 ~ y) (/ T 2ma 2 ) 1 

2m (1+y) J 

Thus, the faster the body is spinning, the faster the yo-yo 
will stop it. The maximum tension in the cord (which is 
assumed massless) is 

r -M-sr 3m V 

— j o(/ r2nui-*)J 

When the stretch yo-yo is used, the effect of errors in 
the body inertia /, initial spin rate <o n , cord length fi and 
despin mass m is minimized. 
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